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Abstract 

The controversial problem of the evolution of an isolated system with an 
internal adiabatic wall is investigated with the use of a simple microscopic 
model and the Boltzmann equation. In the case of two infinite volume one- 
dimensional ideal fluids separated by a piston whose mass is equal to the mass 
of the fluid particles we obtain a rigorous explicit stationary non-equilibrium 
solution of the Boltzmann equation. It is shown that at equal pressures on 
both sides of the piston, the temperature difference induces a non-zero average 
velocity, oriented toward the region of higher temperature. It thus turns out 
that despite the absence of macroscopic forces the asymmetry of fluctuations 
results in a systematic macroscopic motion. This remarkable effect is analo- 
gous to the dynamics of stochastic ratchets, where fluctuations conspire with 
spatial anisotropy to generate directed motion. However, a different mecha- 
nism is involved here. The relevance of the discovered motion to the adiabatic 
piston problem is discussed. 
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I. INTRODUCTION 



The recent work of E.Lieb and J.Yngvason [l|] on the foundations of thermostatics has 
given a new interest to an old but still controversial problem of thermodynamics. This is 
the so called "adiabatic piston problem" described as follows (see also recent references 
in 0). The system is composed of a finite volume cylinder containing two fiuids separated 
by an adiabatic movable piston, as shown in Fig.l. A brake maintains initially the piston at 
rest, and the two fiuids are in equilibrium with temperature, pressure and volume given by 
(Ti(°\pf\ = A(Li + XW)), and (Tf , ^^^^ = ^(Ls-XW)), respectively. A denotes 
here the cross section area of the cylinder, and X^°) the initial position of the piston. At a 
certain time the brake is released, and the problem is to find the final equilibrium state. 
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Figure 1: Adiabatic piston problem 

The laws of thermostatics give immediately the condition for the mechanical equilibrium, 
namely that in the final equilibrium state the pressures must be equal, i.e. p{ = "p^ (with 
E{ + Ei = + E^). However, they give no information concerning the final temperatures, 
except for the condition that they must be such that for both fiuids the final entropy cannot 
be smaller than its initial value, i.e. Sl > 5*°, and > 82- 

On the other hand, from the laws of thermodynamics one can derive ordinary differential 
equations for the time evolution of the system. Solving these equations would then give the 
unique solution of the problem, and in particular the final temperatures. This idea has 
been recently followed in ||^, where a system of thermodynamic evolution equations has 
been derived within a very simple model. One had to distinguish then between the case 
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where the two fluids have vanishing viscosities and the case where at least one of the fluids 
has a non-zero viscosity. For non-viscous fluids the thermodynamic equations predicts a 
periodic motion subject to the energy conservation law, with constant subsystem entropies: 
S^{t) = S2{t) = Si°\ When at least one of the fluids is viscous, the system will 

attain the flnal equilibrium state with equal pressures pi = P2- However, to flnd the flnal 
entropies (or temperatures) one has to solve the evolution equations which so far could not 
be achieved explicitly. These same evolution equations were also obtained in the framework 
of the kinetic theory of gases, assuming the fluctuations of the piston to be negligibly small 
and considering perfect gases, and several numerical solutions have been obtained 0. 

But even having obtained the macroscopic description of the evolution of the system 
the controversial question remains whether the stochastic frictionless motion of the piston 
induced by collisions with the fluid particles could eventually lead (on a much longer time 
scale) to a flnal state where the temperatures of the fluids are equal, and whether such a 
microscopic fluctuation mechanism would not violate the second law of thermodynamics. 
Clearly, the analysis of such problem requires the microscopic description of the dynamical 
evolution. 

In this article we thus present a microscopic study of the motion of the adiabatic piston in 
contact with two inflnite volume reservoirs. Our analysis, based on the Boltzmann equation, 
has been inspired by the above fundamental adiabatic piston question. 

In Section 2 we deflne the microscopic model and we discuss the applicability of the 
Boltzmann equation to the description of the motion of the piston. The rigorous determi- 
nation of the stationary state in one dimension for the case where the mass of the piston 
is equal to the fluid particle mass is presented in Section 3. The derived solution at equal 
temperatures reduces simply to the Maxwell equilibrium state. But when the temperatures 
are different a non-equilibrium stationary state appears characterized by an average velocity 
of the piston oriented toward the higher temperature region, although no macroscopic force 
is present (the pressures on both sides of the piston are supposed to be equal). This situation 
is thus precisely that of the "molecular motor problem". The rigorous derivation (within 
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Boltzmann's theory) of the fluctuation induced macroscopic motion is the main result of our 
work. The relevance of this prediction to a much more realistic problem of a massive piston 
is discussed in the closing Section 4. The case where the mass of the piston is much larger 
than that of a fluid particle will be discussed in the forthcoming paper |Q . 

II. ADIABATIC PISTON SEPARATING TWO THERMOSTATS: THE 

BOLTZMANN EQUATION 

In order to simplify the analysis as much as possible we consider in the following the 
one-dimensional system occupying the whole infinite line. The two fluids are assumed to 
be perfect gases composed of identical point particles of mass m. They are separated by a 
structureless particle of mass M which plays the role of the piston. The perfectly elastic 
collisions between the fluid particles result in instantaneous exchanges of velocities. So, from 
the point of view of the dynamics of the piston the fluid particles can be considered as non- 
interacting, penetrating each other in their unperturbed free motion. The only perturbation 
of the fluid states is due to collisions with the piston. The similar model has been already 
investigated almost forty years ago by Lebowitz and few years later by Alkemade, van 
Kampen and Mac Donald [|1^ . 

When a collision occurs, the pre-collisional velocities V and v of the piston and the fluid 
particle, respectively, suffer a change according to the elastic collision law 

V^V' = V + — {v-V) (1) 

2M 

V V = V (V — V I 

M + m ' 

Let us suppose that at the initial moment the fluids to the left and to the right of the piston 
are at equilibrium with temperatures Ti and T2, and number densities ni and n2, respectively. 
Moreover, the macroscopic mechanical equilibrium will be assumed: the pressures on both 
sides of the piston are equal 

pi = riikBTi = P2= n2kBT2 (2) 



{ks is Boltzmann's constant). 

When M ^ m, the probabihty that the piston interacts back with the fluid particle it 
already collided with becomes negligible. In the absence of recoUisions the piston always 
"sees" the unperturbed Maxwell distributions of the arriving fluid particles, and no correla- 
tions can occur between its velocity and the precoUisional velocity of the fluid particles. The 
perturbations caused by collisions with the piston go away to infinity. So, in the large mass 
limit m/M — > 0, we can expect the linear Boltzmann equation to give a correct description 
of the dynamical evolution of the piston. In order to write down this equation we denote by 
4>^^ and the Maxwell velocity distributions of the two surrounding fiuids, where 

Let F{X, V; t) denote the probability density for finding the piston at point X with velocity 
V at time t. When m/M <^ 1, F{X,V;t) satisfies the linear Boltzmann equation of the 
form 

di + ^dx) ^' ^ / ^''l^ ~ 

-e{v-V)n,<P^^{v)F{X, V;t) 

-e{V -v)n24>t{v)F{X,V-t)] 

where 9 is the unit Heaviside step function. 

Our aim is to determine the velocity distribution of the piston in the asymptotic station- 
ary state attained for i — > oo. It turns out that this problem can be given a rigorous answer 
in the special case of M = m. But when the mass M of the piston is equal to the fluid 
particle mass m, the recoUisions can take place, and the validity of the Boltzmann equation 
becomes questionable. However, what we are really interested in is the case of M » m. We 
thus expect that neglecting the efi'ect of recoUisions in the equal mass case will lead to the 
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evolution qualitatively similar to that occuring for M ^ m, simply because the same class 
of dynamical events will be taken into account. This should apply in particular to the nature 
of the stationary state. It should be noticed here that even when compared with the exact 
results the linear Boltzmann equation for M = m represented a very reasonable approxima- 
tion to the velocity correlation function and to the Van Hove self-correlation function of the 
one-dimensional hard rod fluid 0. 

Having adopted the above point of view we are going to study in the next section the 
equation satisfied by the stationary probability density ^{V) for finding the piston with 
velocity V in the equal pressure case (see eq.(0)). Denoting by F{X,V) the stationary 
solution of equation we have 

^V) = J dXF{X, V) (5) 



III. SOLVING THE BOLTZMANN EQUATION: APPEARANCE OF A 

MOLECULAR MOTOR 



According to the linear Boltzmann equation @ the stationary velocity distribution ^{V) 
in the equal pressure case (||) satisfies the equation 



dv\V — v\ 



-1-1 J-2 



(6) 



dv\V — v\ 



e{y - v)^ct^^^{y) + e{v - v)^ct^^^{y) 

J-l J-2 



In the following we shall assume T2 > Ti. 

It is convenient to express eq.(^ in terms of the following two functions 



m 

G{V) 



V 

00 

V 



dv{v-V)^<f>J?^{v)+ r dv{V-v)^<f>Jp^{v) 
1 1 J —00 1 2 

dv{v - V)^{v) + r dv{V - v)^{v) 



(7) 
(8) 



Using the identities 



\V - v\e{V -v) = ]^[{V -v) + \V -v\ 



\v - v\e{v -v) = -v) + \v -v\ 



we rewrite eq.(^ in the form 
I{VMV) - \[^ct^^{V) + ^r^,{V)]G{V) = < V >)[^<P^,{V) - ^^'^,{V)] (9) 

Z il I2 I ll ±2 

where <V > denotes the average velocity of the piston 

<V>= J dV^{V)V (10) 

We now notice that the first- and second-order derivatives of the functions I{V) and G{V) 
are given by 

I'(V) = - / dv-<l,J?^iv) + / dv-<l>^^iv) (11) 

JV 1 1 J —00 1 2 

I"{V) = ^<P^^{v) + ^<P^^^{v) (12) 

J 1 -1-2 

POO rV 

G'{V) = - dv<!>{v) + / dv<!>{v) (13) 
G"{V) = 2<^{V) (14) 

Therefore takes the form of a simple differential equation 

I{V)G"{V) - I"{V)G{V) = [I{V)G'{V) - I'{V)G{V)]' (15) 



iV-<V>)[^<t>^^iV)-^<l>J^^iV)] 

J-1 J-2 



Using the asymptotic formulae 



. T^V, for V +00 

(16) 

-^V, for V ^-00 



, iV- < V >), for V ^+00 
G{V)^{ (17) 

-{V- < V >), for V ^-00 



we can readily evaluate the limit 



^limJI{V)G'{V) - I'{V)G{V)] = + ^) < ^ > 
Integrating equation (|T5|) with the boundary condition (|18D yields the relation 



(18) 



I{V)G'{V) - I'{V)G{V) = '^mV) - <Pnin 



+ <V> [j^ dv^c^^^iv) + f_Jv^J>^^{v)] 



(19) 



We thus arrived at a first order differential equation for the unknown function G. The 
solution of eq.(|l^) can be now constructed by standard methods. Introducing the auxiliary 
function H{y) 



G{V)=I{V)H{V), 



(20) 



and using the relations 

lim H(V) = lim = T2 

lim H{V)= lim ^^ = Ti 

V^-OD ^ ' V^-oo I(y) 

one readily determines the function B. . It reads 



(21) 



(22) 



+ 



< > 



Here sgn(K — u) = 6{V — u) — 6{u — V). The relations (0), ( PUD permit to determine 
the shape of the stationary probability density ^{V). We find 



$(F) 



{V- <V>] 



J-l J-2 



(23) 



TI + T2 + 



— ^sgn(V-t.) ^(0- («))+< V> / dv-<l>Jp^iv) + dv-<l>^^iv] 



In order to complete the calculation of ^{V) we have still to evaluate the average piston 
velocity < > of the piston which shows in equation (pSf). This can be achieved from the 
normalization condition 



J dV<^{V) = 1 



(24) 



To perform the calculation we group together the terms multiplying < V >, and rewrite 
equation (|23| ) in the form 



K{V) <V>= L{V) 



(25) 



It turns out that the function K{V) can be written as a derivative 



K{V) 



d_ 
dV 



y 1 /-oo 1 



We then deduce from fEUf) the relation 



(26) 



u 1 roc 1 

dv-<l>J^^iv) + dv-<l,JF^iv)) (27) 



foo 1/1 1 \ r+oo fill 

dVK{V) = U^-^) yIp. 



which follows directly from the asymptotic formulae ([TgD . The right hand side of equation 
(p5|) equals 



L(v) = 2$(\/)-i(ri + r2) 



J-l J-2 



where 



(28) 



AiV) 



i{vy 



-oo i 2 



dv-<P^^{v) 
V Ti 



dv 
iHv) 



Taking again into account the asymptotics ( |16D and the normalization (^) we get 



dVL{V) 



(Ti - T2 
2T1T2 



T,-T2 + ^ 
m 



dv 
iHv] 



(29) 



Integrating equation (|25|) over the whole velocity space we thus arrive at the following 
formula for the stationary mean velocity of the piston 



<V> 



(30) 



du 



P{u) 



U 1 POO 1 

dvTfr<l>TM+ dv-ct>^^{v) 

oo J 2 ''U 1\ 



The formula (|30| ) is the main result of our analysis. Whereas the mean velocity vanishes 
for Ti = T2, the piston acquires a macroscopic systematic motion when Ti 7^ T2, and this 
under the equal pressures condition (j^). 

In order to further investigate this phenomenon let us analyze the case of a small tem- 
perature difference. We thus put 

ST 

Ti=T, T2 = T + 5T, with < — < 1 



When 5T ^ 0, the denominator in equation (|30|) takes the asymptotic form 



m 



dwh{w) + 0{5T) 



(31) 



where 



h{w) 



dw\u — w\(l){w) 



(32) 



and d)(w) is the dimensionless Maxwell distribution 



(p{w) = 

On the other hand the numerator equals 

1 



exp(— ti?^/2) 



ST 



1 — - / dwh{w){w'^ — l)(f){w) 



(33) 



So, when ST — » 0, the mean velocity of the piston is given by 



<V>-- 



iknT ST 



m \ T 



1 — - / dwh{w){w'^ — l)0(w) 



dwh{w) 



^ -1 



+ 0{ST) (34) 



Integration by parts yields the relation 
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J dwh{w){w'^ — = J dw w(f){w)h'{w) 

The inequality 

wh'{w) = -2[h{w)f^'^w J dw'sign{w - w')(j){w') < (35) 

permits to determine the sign of the mean velocity of the piston. We find 

< > > 0, if 5T > 0. (36) 

The inequality ( ^61) shows that in the stationary state the piston performs an uniform 
motion in the direction of the higher temperature region. In fact, we checked by numerical 
analysis that the inequality (|36|) remained true for any T2 > Ti. 



IV. CONCLUDING COMMENTS 

We have studied the effect of fluctuations on the motion of the piston separating two 
semi-infinite regions filled with perfect gases at equal pressures. The explicit stationary 
solution of the Boltzmann equation (p3D, ( |30D obtained for the special case where the mass 
of the piston equals that of a fluid particle reduces to the equilibrium Maxwell distribution 
when the temperatures on both sides of the piston are equal. The mean velocity then 
vanishes. In fact, it has been rigorously shown that the piston (tagged particle) performs in 
this clasical diffusive motion [Hi, HI. 

When the temperatures are different, the space asymmetry induces asymmetry of fluc- 
tuations which turns out to have a macroscopic effect. The piston attains a stationary state 
with a non-zero average velocity oriented toward the higher temperature region. In view 
of the absence of macroscopic forces we discover here a genuine molecular motor. There is 
some analogy between our prediction and the "stochastic ratchet" problem where one ob- 
serves that fluctuations conspire with the spatial anisotropy to induce a macroscopic directed 
motion in the absence of macroscopic forces [O, W2 . 
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We expect the same qualitative behavior in the case of a massive piston whose mass is 
much larger than that of a gas molecule. Indeed, when M ^ m the Boltzmann equation, 
which was the basis of our analysis, gives an adequate description of the dynamics because 
the creation of dynamical correlations between the piston and the gas molecules via the 
recoUision mechanism becomes negligible. The persistence of the non- vanishing average 
velocity in the M ^ m regime at equal pressures and different temperatures will be the 
subject of the forthcoming article |^. One can thus expect that the stochastic motion of 
the piston as a whole can lead to an effective energy transfer between the separated fluids 
in the absence of pressure difference. Hence, although the piston has no internal degrees of 
freedom calling it adiabatic becomes a misuse of the term. 

It remains to be seen what is the relevance of our analysis to the original adiabatic 
piston problem where the fluids occupy finite volumes. This is a delicate question indeed, 
because in one dimension even for a very massive piston the collisions with the boundaries 
will inevitably induce recollisions: the piston will interact back with the particles it had 
already met before. So, the Boltzmann description cannot be reliable in the long time 
limit. Indeed, at the time scale characterized by the time necessary to cover the total length 
of the cylinder (Li + L2) with the thermal velocity of the fluids the effects of dynamical 
correlations will start to play an important role. However, for sufficiently large volumes, 
after the mechanical equilibrium has been established, one can expect (on a larger time 
scale) the occurence of a net displacement of the piston, reducing the volume of the the 
higher temperature fluid, as it is the case for the conducting piston. It seems likely that for 
finite volumes the Boltzmann description can work much better in two and three dimensions 
because then the recollision effects due to boundaries become weaker. Anyway, for finite 
volumes one will have to consider the system of kinetic equations coupling the evolution of 
the state of the piston to that of the fluids. 
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